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By calculating the change of entropy, we prove that the first law of black hole thermodynamics 
leads to the tunneling probability of massive particles through the horizon, including the tunneling 
probability of massive charged particles from the Reissner-Nordstrom black hole and the Kerr- 
Newman black hole. Novelly, we find the trajectories of massive particles are close to that of 
massless particles near the horizon, although the trajectories of massive charged particles may be 
affected by electromagnetic forces. We show that Hawking radiation as massive particles tunneling 
does not lead to violation of the weak cosmic-censorship conjecture. 
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INTRODUCTION 



In 1974, Hawking discovered [![, when considering quantum effects, that a black hole can emit thermal radiation, 
which was also confirmed by many other scientists [2|-|4j. Hawking's result does not coexist with energy conservation, 
however, since the spacetime background is fixed in Hawking's calculation Recently, Parikh and Wilczek 
developed the method of Hawking radiation as tunneling, in which the spacetime background was not fixed and the 
reaction was included. This method has been discussed in different situations [7 - 11 1 , and its self-consistency has been 



checked by using thermodynamic relations [12|-|14| 



The laws of thermodynamics play an important role in the area of study on the black hole. Although the method 
l— ~~ '■ of Hawking radiation as tunneling can give the radiation rate, the relation between tunneling dynamics and black 
hole thermodynamics still needs to be checked. Recently, Pilling [p| suggested a method of Hawking radiation as 
tunneling from black hole thermodynamics. In this method, the change of entropy is calculated to obtain the tunneling 
probability directly from the first law of black hole thermodynamics. The connection between black hole tunneling 



0^ ■ and thermodynamics has been verified further 16 1, when the tunneling probability through quantum horizon (a. Ill 
[20j had been obtained directly from the first law of black hole thermodynamics [l6T |. However, these discussions 
only involve the massless particle tunneling. Since the massless particles do not take charge with themselves, the 
trajectory will not be affected by the electromagnetic forces. It is noted that the tunneling probability of massive 
charged particles from the Reissner-Nordstrom black hole 7] or Kerr-Newman black hole [lfj| has been obtained, but 
whether the tunneling probability can be obtained directly from the first law of black hole thermodynamics or not 
• • i has not been considered. On the other hand, since the massive particles can take charge and angular momentum with 
themselves, the weak cosmic-censorship conjecture (21rl23| should be revisited with the tunneling of massive particles. 
In classical theory, the weak cosmic censorship conjecture is able to get full support, but the quantum effect may 
challenge it [2J|. So, it is interesting to verify whether the Hawking radiation as massive particles tunneling could 
lead to violate the weak cosmic-censorship conjecture or not. 

In this paper, we will calculate the massive particles' tunneling probability from the first law of thermodynamics 
and apply this method to the tunneling from Reissner-Nordstrom black hole and Kerr-Newman black hole. Novelly, 
we find the trajectories of massive charged particles are almost identical to those of massless particles near the horizon, 
although the trajectories of the massive charged particles connect with the electromagnetic forces. We will also show 
that the massive particles tunneling of a black hole does not violate the weak cosmic-censorship conjecture, while, 
in some situations, the Reissner-Nordstrom and the Kerr-Newman black hole may evolve into extreme black holes 



The paper is organized as follows. In the second section, we calculate the change of entropy to obtain the probability 
of massive particle tunneling from the first law of black hole thermodynamics for a general class of static, spherically 
symmetric spacetime. In the third section, we investigate the massive charged particle tunneling from the Kerr- 
Newman black hole and then discuss the situation of Reissner-Nordstrom black hole. Next, we discuss the weak 
cosmic censorship conjecture. Finally, we give some discussion and conclusions. 

In this paper we take the unit convention k = h = c = G= 1. 
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MASSIVE PARTICLES' TUNNELING AND THERMODYNAMICS 



In this section, we will extend the method by Pilling [15[ to the situation of the massive particle tunneling from 
black hole thermodynamics. Let us start with the metric 



ds 2 = -f(r)dt 2 s + -— +rW, (1) 

where the metric describes a general class of static, spherically symmetric spacetime. It is noted that there is a 
coordinate singularity in the metric at the horizon r = R which is given by / (R) — g (R) = 0. For the sake of 
convenience, we can make the Painleve coordinate transformation to remove the singularity at the horizon, 



dt^dt s -J ) g(r) dr, (2) 
V f( r )9{r) 

where the time transformation is dependent on r, not t s . So the metric remains stationary and the time direction is 
still a Killing vector. After the transformation, the metric (QJ becomes of the following form: 



ds 2 = -f (r) dt 2 + 2/ (r) J ], f^\ dtdr + dr 2 + rW. (3) 

V f{ r )g( r ) 

The new metric is called a Painleve metric that is regular at the horizon. The time t here is that of radially free-falling 
observer through the horizon. 

Generally, for the static coordinates, the coordinate time can represent the global time due to Einstein's simultaneity. 
For example, the Schwarzschild coordinate time in global space can be regarded as the time recorded by a standard 
clock that is resting at spatial infinity. The coordinates © is stationary and not static (the component goi of the 
metric tensor does not vanish in the metric ©), so we must check whether the time is global. According to the general 
theory [28| of the coordinate clock synchronization in a spacetime, made by Laudau and Lifshitz, the condition that 
the simultaneity of coordinate clocks can be transmitted from one place to another is either the integral / f^^' = 
or the integral J —dx l is independent on the path. The condition can be described as another form Q equivalently, 
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Obviously, the coordinate © satisfies the condition so that the coordinate time in Painleve coordinates can be used 
as the global time and the coordinate clock synchronization can be transmitted from one place to another. This is 
important for tunneling process in quantum mechanics. 

From the coordinate ((3]), the radial null geodesies can be obtained as 




where the positive (negative) sign gives the outgoing (incoming) radial geodesic, under the implicit assumption that 
t increases towards the future. 

When considering the massive test particles moving radially in the background, the radial null geodesies is improper 
for such particles. In analogy to Ref. p}, we treat the massive particle as de Broglie wave and obtain the expression 
for r, 



1 1.9oo 1 / 

r = Vp= 2^ = -2gVr^ 
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where v p is the phase velocity, and v g is the group velocity. It is obvious that the geodesies of massive particles are 
different from radial null geodesies. The surface gravity of the black hole for the transformed metric ([3]) at the horizon 
is given by 



1 l-g(r) ,df(r)\ 

Generally, when discussing the surface gravity of a black hole, we are defining a notion that behaves analogously to 
the Newtonian surface gravity, while these two things are not the same because the acceleration of a test body at the 
event horizon of a black hole turns out to be infinite in relativity Generally, the surface gravity of a black hole is 
not well defined. However, we can define the surface gravity for a black hole whose event horizon is a Killing horizon. 
According to the equation in terms of time-like Killing K as K a S7bK a — kqKi,, we can obtain the surface gravity 
in terms of Christoffel components as (JT])- Note that, in static spacetime, a more general expression for black hole 
surface gravity should be taken as the geometrical surface gravity [29| , 

K = 7K (8) 

which gives an additional factor comparing with the Eq. ([7]) and this factor is a constant at horizon. 
Due to / (R) = g (R) = 0, we can expand / (R) and g (R) near the horizon in powers of r — R, 

f(r) = f (R) (r-R)+o((r-R) 2 y 

g (r) = g (R) (r - R) + O ((r - R) 2 ) . (9) 

According to the laws of black hole thermodynamics, the Hawking temperature is expressed in terms of the surface 
gravity via 



k y/f (R)g' (R) . . 

Th= 2^- 4^ ■ (10) 

Note that the Eq. (flOl) is not general local temperature expression at horizon, since it lacks the relative factor stemmed 
from the use of Kodama vector instead of the static Killing vector [29( . The general expression for the temperature 
can be taken as 

T= — = —. 11) 



Thus we must treat the first law of black hole thermodynamics and the entropy carefully. The entropy should be 

A _ n 

4-y J 



given by the area, A, of event horizon as S = -A — —R 2 . Thus the first law of black hole thermodynamics maintains 



its form TdS = dM. 

Let us consider the black hole thermodynamics in the region near the horizon. When the mass of the black hole 
changes from Mi to Mf, the change of the entropy is given as 

^=[dS=[ Mf ^-dM=[ Mf ^R^dM. (12) 
J J Mi dM J Mi 7 dM V ' 

Considering the small path near R, we can insert the mathematical identity Im -^z^dr — — n in the formula (|12l) . 
Thus we obtain 

2t f M f rrj R dR 



AS = —Im / ———dM. (13) 

7 JMi Jr, r~R dM 



With (jlip and the expression of the temperature in thermodynamics y = we can get 
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R dR = 2 
dM ~ V/' (R)9 W 

Using the Taylor series ©, we can gain the geodesic ^ of a massive particle near the horizon, 



r = -\\[f' ( R )9 (R) (r-R) + o((r- i?) 2 ) . (15) 

It is noticed that the geodesic of massive particle is equivalent to that of massless particle [IH in the first order. That 
is to say, no matter whether the particles have mass or not, their trajectories will be close to each other when they 
approach the horizon. In particular, their geodesies coincide precisely with each other at the horizon. One reason 
for this phenomenon is that the wavelength of the outgoing particle is infinitely blueshifted at the horizon. This is 
also the key argument that the WKB approximation can be used for Hawking radiation as tunneling. Sometimes, 
the infinite blueshift property of the event horizon could lead to the so-called "trans-Planckian puzzle" [3(|, which 
is considered as a mathematical artifact of horizon calculations nowadays [31] . On the other hand, the tunneling 
process is instantaneous, so the outgoing particles, both the massless and massive particles, take the same quickest 
path across the horizon. 

With equations (Qj} and (TT5|) , we can obtain the final form of the change of entropy (fl~3l) as 



2 f M f r r f dR 2 
AS = Im / / — dM = Iml. (16) 

7 JMi Jr, r 7 

where I is the action for an s-wave outgoing positive particle in WKB approximation. Consequently, the tunneling 
probability is given as 



r ~ e~^ Im/ = e AS . (17) 
Thus we obtain the massive particles tunneling probability from the change of entropy as a direct consequence of the 



first law of black hole thermodynamics, along the line presented in Ref. 15j. Note that when 7 = 1, the entropy 



obeys the Bekenstein-Hawking expression. But the black hole entropy may have other expressions, for example, when 
7 = j, one has 5=4, which can be obtained from black hole solution in modified f(R) gravity (32[. This implies 
that it isn't always true in general that the related black hole entropy is given by the Bekenstein-Hawking area law 
32] , while the laws of black hole thermodynamics play an important role in Pilling's argument [TBI j . 

THERMODYNAMICS ABOUT KERR-NEWMAN BLACK HOLE 

The Painleve line element of the Kerr-Newman black hole is 



AS - E , 2 , J{2Mr~Q 2 )(r 2 + a 2 )^ 



ds z = ^ ^-dt z + — -rir^ + 2 v v — ' V dtdr + Ed6 2 , (18) 

(r 2 + a 2 ) 2 - Aa 2 sin 2 6 r 2 + a 2 ( r 2 + a 2 ) 2 - Aa 2 sin 2 6 

Qr (r 2 +a 2 ) 

A = , 2 2^ 2 A 2 2 dt = Adt > (19) 

(r 2 +a 2 ) - Aa 2 sm 2 6 

which is obtained from the Kerr-Newman black hole in the Boyer-Lindquist coordinate system by the generalized 
Painleve-type coordinate transformation, 



^-^-iV^^*- (20) 
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and the dragging coordinate transformation 



a (r 2 + a 2 - A) 
(r 2 + a 2 f - Aa 2 sin 2 9 



# = — 5 ^-<ft, (21) 



where E = r 2 + a 2 cos 2 9, A = r 2 + a 2 + Q 2 — 2Mr, and </> fc is the coordinates before transformation. The form of 
the vector potential A is unchanged up to a gauge transformation. A t is the electric potential, and the parameters 
M, Q and J = Ma are the mass, the electric charge, and the angular momentum of the black hole, respectively. The 
dragging coordinate transformation is made in order to make the event horizon coincide with the infinite redshift 
surface. So the geometrical optical limit can be applied. 

The radial time-like geodesies of massive charged particles are given by 

1 l0« A 

r = v v = ~ v a = ~ ~ — = , (221 

P 2 9 2g tr 2y/(r 2 + a 2 )(r 2 +a 2 -A) 



Near the outer horizon r = i? + = M + y ^ 2 — Q 2 — a 2 , which can be obtained by solving A = 0, and the geodesies 
can be expanded in powers of r — R + as 



(M + y/M 2 -Q 2 - a 2 ) + a 
The radiation temperature is obtained as 



^ M 2 - Q2 - a2 <r -R+) + Ur - R+f) . (23) 



« 1 ^M 2 -Q 2 -a 2 

H= 9^ = 9^7 — ~ • (24) 

' M + y/M 2 -Q 2 -a 2 

Let us consider the black hole thermodynamics in the region near the horizon. If the energy of the black hole 
changes from Hi to Hf, the change of the entropy S = ir (i? +2 + a 2 ) is 

. _ f a * ( dS ,„ r dS „ dS ,\ 

AS = J Hi {dM dM+ dQ dQ+ dJ dJ )- (25) 
It can be verified that Eq. (|25p is another expression of the first law of black hole thermodynamics, dM = T^dS + 

Qr(r 2 +a 2 ) 

<&dQ + ildJ, where $ is the electric potential which is defined by $ = A t = ^ 2+a 2^_ Aa 2 g in 2 g > an d is the dragging 
angular velocity which is defined by ft = ^ = = a ^I,^ +a , , A ~? 7 „ . At the horizon, we have $ + = $L-r+ = 

° J J at (r 2 +a 2 ) -Aa 2 sin 2 ' ir—n^ 

q(m+^J M 2 -Q 2 -aA 

2-^ — , 51 + = 57| r=fl + = -r a — r-2 . Considering the conservation of energy, charge and 



^M+y/M 2 -Q 2 -a 2 ^ +a 2 ' --.<"' ^ M+y/M 2 -Q 2 -a 2 ) +a 2 

angular momentum, we have 



dS 




_dR+ 




= 2ttR^ 




dM 


dM 


dS 


= 2ttR a 


_8R+ 






dQ 


dQ 


dS 


= 2ttR a 
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2tt [M + y/M 2 - Q 2 - a- 
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a 



y/M 2 -Q 2 -a 2 



-2-kQ [M + y/M 2 -Q 2 -d< 

y/M 2 - Q 2 -H? 

— 2-na 



dJ y/M 2 - Q 2 - a 2 ' 



(26) 

(27) 
(28) 
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where a = J/M. Similarly, these expressions (|26[) . (|27|) and (|28)l are consistent with the first law of black hole 
thermodynamics, that is, = ^2L, ^ = and 37 = ^f-^ + - Considering the small path near the horizon, we 

can insert the mathematical identity Im J^f r _ 1 fi+ dr = —tt in the formula ([2"5jl . Thus we obtain 

A5 = Im / / -dr - — + -^dQ + ^—dJ 

* Jm Jn \r~R+ )\dM dQ V dJ 

rHf rrf / , $+dr f2+dr \ 

= _2Im / / -. a dM - \, dQ - \a A > (29) 

Jm Jn \K(r-R+) n(r-R+) * n(r-R+) )> K } 

where we have used the expressions of the temperature in thermodynamics ~ = Jp-, the electric potential <£> + and 
the angular momentum fi+. Then Eqs. (|23|) and ([26)) - (f28| give the form of the change of entropy as 

fHf frf / j j j \ 

AS = -2Im / — dM-$+— dQ-n+— dJ) (30) 
VjJi Jn \r r r ) 

Using the Hamilton's equations, r = Igrl^,^,^) = §r, A = M;\(A*r> p "***) = and ^ = 

4p-\r±, pap ) — ^ + irp~-> we can express the change of entropy in an explicit form as 



pr f 

AS = -21m / (P r r - P At A t - dt = -2ImJ (31) 

J ri 

which is related to the emission rate of tunneling particle by T ~ e~ 2JmI = e AS . 

Thus we obtain the tunneling probability from the change of entropy of Kerr-Newman black hole as a direct 
consequence of the first law of black hole thermodynamics. It is noticed that when considering the conservation of 
energy, charge and angular momentum, the tunneling particle can take the charge and angular momentum. Note that 
the generalized coordinate At is an ignorable one and the coordinate <f> is a cyclic one, so the action can be written as 
fj (p r r-P At A t -P^) dt. 

In Ref. [10], it has been pointed out that the tunneling process is closely related to the first law of black hole 
thermodynamics, but the change of angular momentum was treated by the change of mass with the relation 5 J = adM, 
i.e., the angular momentum of the unit mass is kept as a constant. In our calculation, although the specific angular 
momentum a is not kept as a constant, the first law of black hole thermodynamics can still lead to the tunneling 
probability by calculating the change of entropy. So, our result seems to be more reasonable and more general. 
Another potential problem is that of how the angular momentum of black hole is carried away in the semiclassical 
WKB approximation, where only s-wave outgoing particles were considered. That is, the potential I (I + 1) /r 2 related 
to the orbital angular momentum does not work when calculating the tunneling probability. The angular momentum 
taken by the outgoing particles with themselves only means rotation around the center of the spherically symmetric 
black hole, not including the self-spinning. That is the reason why the specific angular momentum is kept as a constant 
m Ref. [lOt, but 

our calculation shows that this restriction is not necessary. 



WEAK COSMIC CENSORSHIP CONJECTURE 



The uniqueness theorem [2l|, [22[ states that all stationary black hole solutions of Einstein-Maxwell equations are 
uniquely determined by three conserved parameters: the gravitational mass M, the electric charge Q, and the angular 
momentum J. The three parameters must satisfy the relation M 2 ^ Q 2 + (J/M) 2 to maintain the black hole as the 
condition demanded by the weak cosmic-censorship conjecture (23[, which asserts that spacetime singularities coming 
from the completely gravitational collapse of a body must be encompassed by the horizon of a black hole. Particularly, 



M 2 = Q 2 + (J/M) 2 characterizes extreme black holes [25H27] . whereas M 2 < Q 2 + (J/M) 2 are concerned with naked 



singularities rather than black holes. In what follows we will check whether the weak cosmic-censorship conjecture 
will hold or not in the situations of massive particles' tunneling. 

For Kerr-Newman black holes, if the angular momentum of unit mass is considered as a constant in the tunneling 
process, the total angular momentum of black holes is not considered as independent variable and its change is closely 
related to the change of mass. Thus at some time the charge had been taken out of the black hole completely and the 
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mass continues to decrease until the extremal case a 2 = M 2 appears, so the temperature is absolutely zero and the 
radiation vanishes. In this case, the weak cosmic-censorship conjecture will not be violated. If the angular momentum 
of unit mass a is not considered as a constant in the tunneling process, the angular momentum and the charge will 
be taken out of the black hole completely before the black hole vanishes because the particles may take angular 
momentum with themselves together with that of rotating around the center of the spherically symmetric black hole. 
For example, the electron's a and Q (suitably specified in geometrized units) both exceed its mass M. In such a 
situation the relation M 2 > Q 2 + a 2 holds until the black hole vanishes and the extremal black hole will not appear. 

Specially, when a = 0, the coordinate (fT8|) decays into the Painleve line element of the Reissner-Nordstrom black 
hole. One can easily approve that the first law of black hole thermodynamics can still give the massive particles' 
tunneling rate. The consideration of massive particles' tunneling is necessary since the massless particles' tunneling 
may lead to the violation of weak cosmic-censorship conjecture. On one hand, the massless particles do not take 
the charges with themselves but take the energy by the self-gravity effect when tunneling out of the black hole. The 
extreme case Q 2 = M 2 can always be reached when the radiation temperature vanishes. On the other hand, there 
exist some massless neutral scalar particles that are spinning. When they tunnel outward, the black hole will rotate 
towards the inverse direction due to the conservation of angular momentum. In this case, no matter how small the 
angular momentum is, the relation M 2 < Q 2 + (J/M) 2 can always be reached, because the particles do not take the 
charge with themselves. Thus, the weak cosmic-censorship conjecture will be violated. It is similar to the case of 
overspinning a nearly extreme charged black hole via a quantum tunneling process [24j . The violation of weak cosmic 
censorship conjecture can be avoided by considering the massive charged particles tunneling. Generally speaking, the 
particle charge-mass ratio may be bigger than one; for example, for the electron e/m ~ 10 11 3> 1, for the proton 
e/m ~ 10 8 3> 1. Thus, the charge will be carried completely out of black hole soon. The situation of violation of 
weak cosmic censorship conjecture will not appear. 

In general, the tunneling particles may contain the massless, massive, massive charged particles in the tunneling 
process. If the rate of mass loss is quicker than that of charge and specific angular momentum loss, the condition 
M 2 = Q 2 + a 2 may be approached and the extreme black holes appear as the terminal of the tunneling. If the charge 
and the angular moment are always taken out of black holes, the extreme situation will not appear in the tunneling 
process. Fortunately, in any case, the weak cosmic-censorship conjecture will not be violated. 



CONCLUSION 



In summary, the self-consistency of Hawking radiation as massive particles' tunneling has been verified by the laws 
of black hole thermodynamics in our paper. We have showed that the probability of massive particles' tunneling can 
be obtained from the first law of thermodynamics by calculating the change of entropy, and this method has also 
been used to gain the tunneling probability of massive charged particles. We have also obtained the tunneling rate 
for Kerr- Newman black hole directly from the first law of black hole thermodynamics when the angular momentum is 
considered as an independent variable, i.e., the angular momentum of unit mass a is not constant in our calculation. 
Moreover, we have proved that the massive particles will be along the same trajectory as that of massless particles 
when tunneling across the horizon. Finally, we have showed the massive particles' tunneling does not violate the weak 
cosmic censorship, while, in some situations, both a Reissner-Nordstrom black hole and a Kerr-Newman black hole 
can evolve into extreme black holes. 
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